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Abstract 

We propose several procedures for creating new families of integer se- 
quences based on the method of Cantor diagonalization. Then we modify 
and generalize this method. The paper includes explicit formulas for most 
proposed families of integer sequences. 
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1 Introduction 

There are well-known transformations of the integer sequences [1],[2],[3]. Let A 
be the set of integer sequences. We construct an array f2 : {w(i,j)}, in which 
each column is an element of A. We define the sequence w(n) G A as the list of 
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the elements of the array in the Cantor anti-diagonal order [4] . All sequences 
in the paper are numbered starting from 1. New sequence starts with: 
w(l)=«;(l,l) ) 

w(2)=«;(l,2) ) w(3) = «;(2,l) ) 

w(4) =u;(l,3), w(5) = w{2,2), w(6) = w(3, 1), 

The sequence w is a triangular array read by rows: row number k contains k 
elements. From [4], [5] the formula for calculating n is: 

n=^ +j -^ i+j - 2 Ki (1) 

The property of (1) is that rearranging i, j we get the list of elements of the 
array f2 in the Cantor diagonal order. 

The formula (1) is obtained by simple transformation of the classical formula: 

(x + y) 2 + 3x + y 
Z = 2 (2) 

The characteristic of (2) is that the enumeration starts from pair x — 0, y = 
to z = 0. 



Using C. Kimberling formula from A002260, we get inverse functions for (1): 

t(t+l) t 2 +3t + 4 
1 = 71 2 ' 3 = 2 

where t = [ — - - — -J . 

Throughout the paper, sequences will be referred to by their Annnnnn num- 
bers, as found in the On-Line Encyclopedia of Integer Sequences [6] . 

We denote the set of integers numbers Z, the set of the integer nonnegative 
numbers Z*, and the set of the integer positive numbers Z+. 



2 Transformation of a Single Sequence 

2.1 Replication 
Example 2.1.1 

Let a £ A: ai,a 2 ,a 3 , ... We form the array f2 replicating the column a. It 
follows that w(i,j) — at for any j: 

ai ai ai 

a 2 a 2 a 2 . . . 

a 3 a 3 a 3 . . . ( 3 ) 
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and u)(n) is formed as: 

Oi, 

ai,a 2 , 

0. 1.0.2, «3, 

The resulting w is a fractal sequence [7] . 

Definition 1. A sequence (3 is called a reluctant sequence of sequence a, if 
(3 is triangle array read by rows: row number k coincides with first k elements 
of the sequence a. 

Some of reluctant sequences are presented in OEIS A002260, A002262, A037126, 
A059268, A215026. 

Formula for a reluctant sequence is: 

Lo(n) = a m , m = n , where t = [ - J . 

Example 2.1.2 

We apply the procedure of replication to the sequence to from Example 2.1.1. 
Using double replication, we obtain to , wich is a reluctant sequence of the se- 
quence lj : 

ai, 

ai,oi, 
a\, oi, a2, 
ai, oi, a 2 , oi, 

01, Oi, a2, Oi, a2, 
Oi, Oi, a2, Oi, a2, 03, 

, , t{t +1) 
cj(n) = a m , m = n , 

where 

t(t + l) , . V8n' — 7 — 1 , . \/8n — 7 — 1 , 
" =« g — ' t=L ^ J, <= L ^ J- 

Example 2.1.3 

We form the array replicating the row a, so w(i,j) — aj for any i. The array 
uj(n) is the transposition of the array from Example 2.1.1: 

ai, 

a 2 ,ai, 

Ct3, 02, 0,1, 



Definition 2. A sequence f3 is called a reverse reluctant sequence of sequence 
a, if (3 is a triangle array read by rows: row number k lists first A: elements of 
the sequence a in reverse order. 
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There are some examples of reverse reluctant sequences - A004736, A025581, 
A104762, A104887, A130321. 

Formula for a reverse reluctant sequence is: 

, , t 2 + 3t + 4 V8n - 7 - 1 
w(n) = a m , to = n, where t = [ J . 

2.2 Self-composition 

Let a £ A: ai,a 2 ,a 3 ,... We will assume in this section a(n) e Z+ for any 
n. The first column of the array is the sequence a: a(n), the second col- 
umn is a („), the third column is a 0o(n) , and so on. Using another notation, 
w(M') = «(«(• ••«(*)))• 

s v ' 

j 

We get new sequence ui(n): 
o>i, a ai ,a2, a aai ,a a2 ,a 3 , 

Example 2.2.1 

Let to e Z+, m > 1 and a(n) = m™, 
New sequence w(n) is: 

TO 1 , TO™ 1 , TO 2 , TO™" ,TO m2 ,TO 3 , 

Example 2.2.2 

Suppose a(n) is the sequence of prime numbers, then the first column of Q is 
the sequence A000040, the second colum of fi is the sequence A006450. The 
start of the sequence w(n) is: 
2, 3, 3, 11, 5, 5, 5381, 31, 11, 7, ... 

Example 2.2.3 

Let <p(n) be Euler's totient function. The first four columns of are sequences 
A000010, A010554, A049099, A049100. The start of the sequence w(n) is: 
1, 1,1, 1,1,2, 1,1,1,2, 1,1,1,1,4, 1,1,1,1,2,2,... 

2.3 Function Acting on an Index 

Let sequence a e A and function /: Z + x Z + — > Z + . We create the array f2: 
j) = a /(iJ) . 

Example 2.3.1 

Let the first column of the the array f2 be the sequence a: a\, ei2, as, ... 

Every next column is formed from previous shifted by fc € Z* elements. We get 

f{hj) = i + kj — k. The start of the sequence w is: 

Gij ak+i,a 2 , a 2 fe+i, a fe+2 , a 3 , a 3fe+ i, a 2 fe+2, a fc+3 , a 4 , 

The sequence is the generalization of A002260 (fc=0), A002024 (fc=l), A128076 
(fc=2), A131914 (fc=3). 
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uj(n) — a m , where m = k(t + 1) + (k — 1) ^ — ^ ^ — nj , 
t-L 2 J- 

Using rearranging i and j in the formula f(i,j) = ki + j — k, we obtain the 
array 0, in which the first row is a, every next row is formed from previous 
shifted by k elements. 

Example 2.3.2 

Let a G A: a\, a 2 , a 3 , k G Z+, and f(i,j)= max {fci + j — k;i + kj — k}. 
The start of the sequence w is: 

Oil a fc+lj a fc+l> a 2fc+l> a fc+2> a 2fc+lj a 3fc+l> a 2fc+2j a 2fc+2i a 3fc+l--- 

The sequence is the generalization of A002024 (fc=l), A204004 (fc=2), A204008 

(fc=3). 

cj(n) = a m , 

where m = k(t + 1) + (k — l)max { - ^ — - — n;n — ^ — }, 

m _ , V8^7-1 , 
t-L 2 J - 

Example 2.3.3 

Let a G .4: Oi, a2, 03, and fc G Z + . The first column of the array f2 is 
the sequence a, in all columns with number j (j > f ) the segment of the se- 
quence a with the length of j — 1: < ak+j-2, flfe+j-3, — , o-k > shifts the sequence 
01, 02, 03, .... 
Therefore 

i — j + 1, if i > j i, 
j — i + k — 1 , if i < j . 



f(i,j) = 



The start of the sequence u> is: 

Olj afc,fl2, Ofe+i, 01,03, Ofc + 2, Ofc, 02, 04, Ofe+3, Ofc+1, Ol, 03 5 05i 

The sequence is the generalization of A143182 (k=2). 
u(n) = a m , 

m ^2 „ , ,,t 2 + 3t + 2-2n, 
where m = |(i + 1) - 2n| + k[ J, 



t+l 



3 Transformation of a Pair of Sequences 

We form the array f2 from two sequences a, f3 G A, a : 01, 02, 03, 

/3 : 61, 62, 63, ••• Let function F: Z x Z — >• Z. We define the array £1 as 
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w (hj) = F(ai, bj). The start of the sequence u> is: 

F(ai,6i), F{a 1 ,b 2 ),F(a 2 ,b 1 ), F(a u 63), F(a 2 ,b 2 ),F{a 3 , h), ... 

Using formulas for reluctant sequence and reverse reluctant sequence we get: 

, t(t + l) t 2 + 3t + 4 
w(n) = F(a mi ,b m2 ), m 1 =n , m 2 = n, 

where t — [ — J . 

We denote to = C{a,0,F). 

Example 3.1 

Let w(i,j) — ai* bj, p and q are primes, a is the sequence p x ,p 2 ,p 3 , 
(3 is the sequence q 1 ,q 2 ,q 3 , ... so u>(n) are numbers of the form p % qK The se- 
quence is the development of A036561. 



Example 3.2 

Let teZ + ,«)(i,j)=a- hence uj(n) are numbers that can be represented 

as the sum of the A:-th powers of elements of sequences a and /3. 

Example 3.3 

We apply the transformation of a pair of sequences for enumeration n-tuple [8] . 
We denote concatenation of two integer positive numbers as 4f . For example 
12-H-345 = 12345. Function F is concatenation: w(i, j) = ai-\\-bj. 
Let r) 1 be 1,2,3,... We denote r] d (n) = C(j? d_1 , t? 1 ,*). 
The array ft for d — 2 is: 

11 12 13 . . . 
21 22 23 . . . 

31 32 33 . . . ( 4 ) 



Starts of sequences r) d (n) are: 

d = 2 11, 12,21, 13,22,31, 14,23,32,41,... (A066686) 

d = 3 111, 112,121, 113,122,211, 114,123,212,131,... 

d = 4 1111, 1112,1121, 1113,1122,1211, 1114,1123,1212,1131,... 

d = h 11111, 11112,11121, 11113,11122,11211, 11114,11123,11212,11131,... 



Example 3.4 

Composition of a pair of sequences. Let ai, bi € Z+ for all i. 
F{a t ,b j ) = b{b{...a i )). 

" v ' 

j 

Start of the sequence u>(n) is: b ai , bb ai , b a2 , bb ba ^ , bb a2 , b a 
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4 Transformation of Several Sequences 



We form the array f2 from I (I G Z + , I > 1) sequences a 1 , a 2 , ...a 1 G A. Elements 
of sequence a 3 denote a{, a 2 , a 3 3 ... Let functions f,s: Z + x Z + — > Z + , 1 < s < I. 
We construct the array 0: w(i,j) = afffjy 

Example 4.1 

Let a 1 , a 2 , ... a' G A The array f2 is the replication of Z columns a 1 , a 2 , ... a 1 : 
f = i, s(i,j) = l+Mod(j - 1,1). 
The start of the sequence u for Z = 3 is: 

a\, al,a,2, a\,a 2 ,a\, a\,a,2,a 2 ,a\, a 2 , a\, a%, a, 2 , a\, 

Using the formula for the reluctant sequence, we get the number of the element 

to of the sequence uj(n) = a r m : 

t(t + l) VF7-1 
to = n , where t = [ J . 

Taking into account the formula from Example 2.1.3, we obtain the number 
of sequence r: 

t 2 4- 3t 4- 4 
r = l+Mod( 3 - n - 1, J), 

where t = I — - I . 

L 2 J 

Example 4.2 

Let a 1 , a 2 , ... a' G A. The array f2 is replication of "a braid" from I sequences: 
f(i,j) = i and = l+Mod(i+j-2,/). 

The start of the sequence w for I = 3 is: 

a 1; a 2 , a 2 , a^aj^a^, a}, a3,, 03, 04, a 2 ,a2,a 2 ,a,4,a 2 , 

Using the formula for the reluctant sequence, we get the number of the ele- 
ment to of the sequence w{n) — a r m : 

i(t+l) V8^7-l 

m = n , where t = . 

2 L 2 J 

Taking into account the formula from Example 2.3.2 at k = 1, we get the 

number of sequence r: 

, , , ,/ ,x , 1 V8n. — 7 — 1 
r = l+Mod(f,Z), where t = [- J. 

Example 4.3 

Let a 1 , a 2 , ... a 1 G A and j G Z + , j > 1. The segment of the sequence a': 
< g4_ 2 , •■•> a i > shifts in all columns with number j the replications of 

sequences a 1 , a 2 , ... a 1 ^ 1 . 
Thereby 

j - », if * < j- 



/(i,j) = 
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f(i,j) is the special case Example 2.3.3 at k = 1. 

Jl+Mod(j-l,/-l), if i>j, 
[/, if«<j- 
The start of the sequence w for Z = 3 is 

a li a lj fl 2i a 2) a l) a 3' a 3' a lJ a 2! a 4i a 4i a 2i a l: a 3: a 5: a 5 j a 3 j a l j a 2 ) a 4 ; a 6 1 ■ • ■ 

Taking into account the formula from Example 2.3.3 at fc = 1, we obtain the 

number of the element m of the sequence cu(n) = a r m : 

,n9 „ , ,t 2 + 3t + 2-2n, , , y/8n - 7 - 1 
m = |(t + l) 2 -2n| + L — J, where * = [- ^ J. 

Formula for number of sequence r is: 

r = l + «(Mod( - + ^ + 4 - n - 1, i - 1) - / + l), 

where i = l ^ 8 "" 7 " 1 ] , v = [(2n - t(t + 1) + l)/(i + 3)J . 

5 Other Pairing Functions 

So far we were using the classical Cantor anti-diagonal order. Next, we define 
several pairing functions Z+ x Z+ — > Z+. 

5.1 Diagonal Enumeration 

We denote < ; W2,(j-i) ; ■•• > a diagonal of the array fi. Any per- 

mutation elements of the diagonal is the enumeration of the array £1 

Example 5.1.1 

We form the permutation so that neighboring diagonals traverse in opposite 
directions. Such sequences were presented previously without explicit formulas 
[9], A056023, A056011. 

Obviously i + j is constant for all elements of the diagonal. We do not change 
the order of the list for an odd i+j. We reverse the diagonal for an even i + j: 

n = \ ((i + j - l)(i + j - 2) - ((-1)*+'' - l)i + + l)j) . 

The result of the enumeration of the array (4) is: 

11, 12,21, 31,22,13, 14,23,32,41, 51,42,33,24,15, 16,25,34,43,52,61,... 
Example 5.1.2 

7 + 1 

We denote by k = l'—^ — J ■ We form the permutation as a symmetrical move- 
ment from center to edges: 

< w k , k ;w(fc-i),(k+i) ;«>(fc+i),(fc-i) ; ••• wi,j ; wj,i > for odd j, 
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< t"(k-i),(fc+i) ; ; ••• «Ji,j ;«Jj,i > for even j. 



i(i + l) + (j-l)(2i+j-4) 

2 5 11 4 — Jl 

^ + l) + 0--l)(2. + i-4) +2( ._. ) _ i; . H< ^ 



2 

The result of the enumeration of the array (4) is: 

11, 12,21, 22,13,31, 23,32,14,41, 33,24,42,15,51, 34,43,25,52,16,61,... 



Example 5.1.3 

7 + 1 

We denote by k = l'—^ — J ■ We form the permutation as a symmetrical move- 
ment from edges to center: 

< wi,j ;wj,i ; W2,(j-i) ; W{j-i),2 ; ••• w k , k > for odd j, 

< w 1 , j ;w j , 1 \ w 2 ,(j-i) ; w (i _i),2 ; ■•• W(k-i),(k+i) \ w (k+i),(k-i) > for even j. 
Such sequences were presented previously without explicit formulas A064578, 
A194982. 

i(2i - 1) + i(j - i)(3* + j - 3), if * < j, 

j(2j - 1) + ^(i - j)(3j + i - 3) + 1, if i > j. 
The result of the enumeration of the array (4) is: 

11, 12,21, 13,31,22, 14,41,23,32, ' 15, 51, 24, 42, 33, 16,61,25,52,34,43, 



Example 5.1.4 

We form the permutation as an asymmetrical movement from edges to center 
and inversely. Neighboring diagonals traverse in opposite directions: 

< wj,i ;w2,(j-i) ;w(j_ 2 ),3 ; ■•• W{j-\),2 ;wi,j > for odd j, 

< ; ,2 ; u>3,(j-2) 5 ■•• W2,(j-i) ;wj,i > for even j. 

n = \ {[i+3 - 1)0 +3 - 2) + ((-l) max { 4 y> + l)i - ((_l)max{i;j} _ X ) A_ 
The result of the enumeration of the array (4) is: 

11, 12,21, 31,22,13, 14,32,23,41, 51,24,33,42,15, 16,52,34,43,25,61,... 
5.2 Angle Enumeration 

We consider squares of the array ft with vertices w(l, 1), to(l, i), w(i, i), w{i, 1). 
We list the array according to the sides of squares from w(l, i) to w(i, i), then 
from w(i,i) to w(i, 1). The numeration was presented previously as A060734, 
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A060736: 

{i 2 — j + 1, if i > j, 
(j-l) 2 + i, i£i<j. 

The sequence w is a triangle array read by rows: row number k contains 2k — 1 
elements. 

Using the angle numeration to the array (4), we have: 

11, 12,22,21, 13,23,33,32,31, 14,24,34,44,43,42,41, ... 

The inverse functions are: 

i = min{t;n- (t - l) 2 }, j = min{t; t 2 - n + 1}, (5) 
where t = [y/n — 1\ + 1. 

Example 5.2.1 

Using the angle numeration to the array from Example 2.3.1, we get: 
ai, afe+i,afe+2,02, a 2 k+i, a,2k+2, a2fe+3, afe +3 , a 3 , 

«3fe+l, 03fe+2, a3fe+3, a3fc+4, «2fe+4, Ofe + 4, a 4 , 

a;(n) = a m , where 

m = (k + i)([V^ TJ -fc^) + t + i, t = n- LV^ - TJ 2 - LV^ - TJ -1. 

Example 5.2.2 

Using the angle numeration to the array from Example 2.3.2, we have: 

O-l, O-k+lj O-k+2, dk+lj fl2fe+l, a2fc+2, 02fe+3j a 2fe+2, a2fc+l, 
a 3k+l, «3fe+2, «3fe+3, 0-3k+4, «3fe+3 , a 3fc+2 , «3fe+l , 

w(n) = a m , where m = (k + 1) [V n — lj + 1 — |n — IV™ — 1J 2 — W n ~ 1J — 1|- 
Example 5.2.3 

Using the angle numeration to the array from Example 2.3.3, we obtain: 

o-i, a / t,a 1 ,a 2 , ak+i,ak,ai 1 a 2l a 3 , a,k+2i «fc+i> o>k, ffli, ct2 ; &3i a 4> 

The sequence is the generalization of A004739 (fc=l), A004738 (fc=2). 

w(n) = a m , where m — kv + (2v — l)(t 2 — n)+t, 
_ n — 1 

t= Lv^ + i/2j + i,«= L— 7-J + 

Using Antonio G. Astudillo formula from A079813, we have formula for v. 
Example 5.2.4 

Stein proposed the boustrophedonic ("ox-plowing") option, "although without 
giving explicit formula" [5]. Natural numbers are placed "by taking clockwise 
and counter-clockwise turns" of A081344. The formula is below: 
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n = 



(i - l) 2 + j + (-ly- 1 ^ - i), if i< j, 
(*- I) 2 + i-(-iy- 1 (i -j), ifi>j. 



The result of the enumeration the array (4) is: 

11, 12,22,21, 31,32,33,23,13, 14,24,34,44,43,42,41, 

51,52,53,54,55,45,35,25,15, ... 



The property of "ox-plowing" is: if t = [V n — lj + 1 is even, i,j are com- 
puted by (5), otherwise i,j are reversed. Therefore we have formulas for the 
inverse functions: 

i=Mod(t; 2)min{t; t 2 - n + 1} + Mod(i + 1; 2)min{i; n - (t - l) 2 }. 
j=Mod(t; 2)min{£; n-(t- l) 2 }+Mod(* + 1; 2)min{f; t 2 - n + 1}, 



6 Generalization of Cantor diagonalization 

The Cantor diagonalization assumes numeration by cells lxl. We consider the 
cover of the array fl by rectangles variable length l m and height h m ; l m ,h m G 
Z+. The upper right corner is covered by rectangle l\ * hi, right to it the rect- 
angle l 2 * hi adjoins, below it the rectangle h * hi adjoins, and so on. Rectangles 
with the length l s are placed in the column, rectangles with the height h r are 
placed in the row. Numeration within all rectangles is row by row, starting from 
first. The rectangle containing the number w(i,j) is called current. 

r s 

We define i?=max{r : h m < i} and S^maxjs : l m < j}. 

m— 1 m—1 

R,Se X* . The length and the height of the current rectangle are ls+i, hn+i, 
respectively. 

Parameter n contains four components: 

• the number of cells in the stair-step triangle with cathetus (R + S) is 

R+S R+S+l-r 

E M E h). 

r—l s=l 

• the number of cells inside R rectangles, lying above the diagonal passing 

R 

through the current rectangle is E W.R+S+2-r- 

r=l 

R 

• the number of the row inside current rectangle before w(i, j) is i— E h r —l. 

r=l 
S 

• the number of the column w(i,j) inside current rectangle is j — E 's- 

8=1 

Therefore 

R+S R+S+l-r R R S 

n =^2 h r( l s)+^2 h rlR+S+2-r + ls+l(-i-^2K-^)+j-^2ls (6) 

r—l s— 1 r—l r—l s— 1 
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If the numeration within rectangles is column by column starting from first, 
formula (6) takes the following form: 

R+S R+S+l-r R S R 

n = E M E l s) + E KlR+S+2-r + h R+1 (j - E Is - 1) + « - E • 
r— 1 s— 1 r— 1 s— 1 r— 1 

Using (2) 

4 _ (i? + ff) 2 + 3i? + 5 
2 

we get the number of the rectangle. Thus we obtain the formula for numer- 
ation within the rectangles row by row or column by column, depending on the 
parity of t: 

R+S R+S+l-r R 

n = E M E U + E KIr+s+2-t + 

r—1 s—1 r—1 

1 ; (ls+i(i - E ^ - 1) + 3 - E i.) - 

f-l)*-l/ s fl \ 

—4: (h R+ i(j- El.-l) + i-EU 

^ v s=l r=l 7 

Suppose t = i? + S 1 ; then the above formula is the numeration within the rectan- 
gles, depending on the parity of the number of diagonal, formed by the rectangles 



Example 6.1 

We apply the enumeration to the array (4) for rectangles with the lengths of 

I, 2,3,... and the heights of 1,2,3,... Numeration inside rectangles is row by row. 
The result is: 

II, 12,13, 21,31, 14,15,16, 22,23,32,33, 41,51,61, ... 

Let the sequences l m , h m be constants: l m = l,h m = h for all m. The example 
I = 2, h = 2 was presented previously without explicit formulas [9]. Using (6), 
we obviously get: 

n = y i( R + S ) 2 + 3R + S) + l(i - hR - S - 1) + j. 

Formulas for R and S take the following form: 
i?=L^Jand5^L 1 y 1 J- 

Example 6.2 

We apply the enumeration for rectangles with the length of 3 and the height of 
2, and the numeration row by row inside it to the array (4). The result is: 
11,12,13,21,22,23 14,15,16,24,25,26, 31,32,33,41,42,43 ... 



12 



7 Transformation of Sequences as Superposition 
of Pairing Functions 

Let /, g be pairing functions Z + x Z + — > Z + , f ^ g. We denote the inverse func- 
tions to / by (fi ,fz ). We create the array w(i,j) = a^-i^ /-!(„)• Using 
function g, we convert the array to the sequence u> : tj(n) = a s (/ _1 (n) / _I (n))" 

Example 7.1 

Using the enumeration from Example 6.2, we place the sequence a, then we 
form the new sequence by enumeration from Example 5.1.2. The start of the 
sequence u> is: 

ax, 02,^4, 05, a 3,ai3, 06, Oi4, 07, Oi6, 015,010,017,08,032, 
026, Oi8, OlOi O33, Og, 035, 
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